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ABSTRACT 

We study the differences in chromospheric structure induced in K stars by steUar 
activity, to expand our previous work for G stars, including the Sun as a star. We 
selected six stars of spectral type K with 0.82< B — V <0.90, including the widely 
studied Epsilon Eridani and a variety of magnetic activity levels. We computed chro- 
mospheric models for the stars in the sample, in most cases in two different moments 
of activity. The models were constructed to obtain the best possible match with the Ca 
II K and the H/3 observed profiles. We also computed in detail the net radiative losses 
for each model to constrain the heating mechanism that can maintain the structure 
in the atmosphere. We find a strong correlation between these losses and Scaii, the 
index generally used as a proxy for activity, as we found for G stars. 
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1 INTRODUCTION 

Solar and stellar chromospheric models have been developed 
to study the dependency of chromospheric plasma parame- 
ters with height and temperature. The best known exam- 
ples are the models for the solar atmosphere com- 
puted by E. Avrett and his co-worker s, in particular 
model C for the average quit Sun by Vernazza et al.l 
(|l98lj) , which was later modified by Fontenla et al. 
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In several cases, these models were used to character- 
ize changes due to activity and spectral type. For example, 
Kelch et al. (1979) studied a sample of eight main-sequence 
stars ranging in spectral type from FO to MO, some of which 
were of similar spectral type and different levels of chromo- 
spheric activity. They computed the photospheric structure 
starting from a radiative equilibrium model for the Teff of 
each star and fitting the Ca 11 K fine wings. The chromo- 
sphere was built using the emission core of the Ca 11 K line. 
To estimate t he radiative cooling rat e in the K line they used 
the Ki index (|Linskv fc Avresll 19781 ) ■ which is calculated as 
the difference between the integrated flux inside the two Ki 
minima of the Ca 11 K line and the corresponding flux for 
the model in radiative equilibrium. 

Their results showed that non-radiative heating is im- 
portant in the lower photosphere of all the late-type stars 
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under study. They found that the value of the Ki index 
and the temperature gradient in the lower chromosphere of 
these stars, as a function of TcH, divides active and inac- 
tive stars, and that the cooling rate in chromospheric lines 
decreases with T^g. Regarding the chromospheric structure, 
they found that the temperature minimum moves outward, 
to lower values of column mass density, with decreasing mag- 
netic activity, i.e. with decreasing non radiative heating in 
the lower chromosphere. 

Semi-empirical models of the dM star AD Leo in its qui- 
escent state and during a flare were built bv lMauas fc Falchil 
(1994, 1996). Subsequently, models of two "basal" (i.e. inac- 
tive) stars of t he same spectral typ e, G1588 and G1628, were 
constructed bv lMauas et al.l ( 19971 ). 

In a previous paper (jVievtes et al.l 2005. hereafter Paper 
1), we computed chromospheric models for a sample of dwarf 
stars of spectral ty pe G, including the Sun as a star, using 
the FAL models iFontenla et al.l (|l993l ) as a starting 
point. Our purpose was to study the changes in chromo- 
spheric structure induced by magnetic activity. The stars 
we modeled were chosen to have similar colors than the Sun, 
and therefore similar photospheric structures, but different 
chromospheric activity levels, probably due to different ages 
and/or rotation periods. These stars can be considered as 
solar analogues, since they share several characteristics with 
the Sun. 

To extend our research to cooler stars and to study how 
the chromospheric structure changes with spectral type and 
chromospheric activity, in this paper we perform a study 
similar to the one in Paper 1 for several dwarfs of spectral 
type K, selected with similar colour, i.e. similar photospheric 
structure, and with different levels of magnetic activity. 
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As the base for our sample we selected one of the most 
studied K stars, Epsilon Eridani (HD 22049), which is an 
active star of spectral type K2 V (-B — V^O.SS), with Tcff = 
5110 K. This star has been widely studied because it is one 
of the ten nearest stars. It has two planets and a belt of 
dust particles around it, which has been compared to the 
Kuiper belt in the Sun. These facts make this stellar system 
resemble our own Solar System. 

Seve ral chromosp heric models have been computed for 
this star. iKelchI (llOTSl ) modeled the lower chromosphere to 
match the Ca II K line profile and integrated fluxes of the 
Mg II h and k lines. Using observations of the ultraviolet 
lines of C II Mg II , Si II and Si III from the lUE satellite, 
ISimon et al.l l|l980l ) obtained a model for Epsilon Eridani, 
which also reproduces hydrogen line profiles not fitted by 
Kelch's model. The thermal structure of this model has the 
onset of the transition zone deeper in the chromosphere and 
a lower temperature in the plateau than Kelch's model. 

A nother chromospheric m odel for Epsilon Eridani is the 
one bv [Thatcher et al.l (llOOll ). who fitted the Ca II K line, 
the infrared triplet line s of Ca II, the Na D doublet, Hq 
and H/3. More recently, ISim fc JordanI (|2005l ). using ultra- 
violet observations from STIS and FUSE, developed a new 
semiempirical model for the upper chromosphere and lower 
transition region of th is star keeping t he ph otosphere and 
lower chromosphere of 'Thatcher et ahl l|l99ll ). 

Finally, Ness fc Jor dan (,20oi )" studied the relative ele- 
ment abundances from the conora and upper transition re- 
gion of Epsilon Eridani, using observations from Chandra, 
EUE, FUSE and XMM-Newton. 

This Paper is arranged as follows: we present our stellar 
sample and discuss the observational data in §2. In §3 we 
describe the models and show the results. In §4 we compute 
the energy requirements to sustain the chromosphere, and 
compare the results with those obtained for G stars in Paper 
I. Finally, in §5 we discuss the results. 



2 OUR STELLAR SAMPLE 

The largest observational study of chromospheric activity is 
the one started in 1966 at the Mount Wilson Observatory, 
which at present includes more than 2200 stars in the spec- 
tral range between F and early K. As the indicator of chro- 
mospheric stellar activity, they use the Scaii index, which 
is the ratio of the fluxes in the H an d K line cores and tw o 
nearby reference windows 20 A wide (|Vaughan et al.iri978l ) . 
The emission in the cores of these lines increase with in- 
creasing chromospheric activity, i.e. with increasing surface 
magnetism. In this work we used the same activity indicator. 

To select the stars in our sample, we require that 
0.82 < B — V < 0.90, a colour similar to e Eri, and that 
the magnetic activity levels are different. All the stars are 
part of the library of southe rn late-type dwarfs published by 
ICincunegui fc Mauad l|2004 hereafter CM04). 

The stellar parameters of the stars in our sample are 
listed in Table [21 In the third column we list the spectral 
type, in the fourth to sixth columns we indicate the colour 
index B — V, T^s and the metallicity. In column 7 we show 
the mean values of the Scaii inde x obtained at the Cerro 
Tololo Inter American Observatory l|Henrv et al.lll99^ , and 
in columns 8 and 9 the maximum and minimum Scaii ob- 
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Figure 1. Scall for each observation for e Eri in our library. The 
open triangles are the values for the different observations, the 
squares indicate the annual averages, and the largest full triangles 
show the two spectra modeled in this paper. 



tained from our spectra (see Cincunegui et al. 2007 for de- 
tails on how this index is obtained) and from the models 
we built in this paper. Finally, in the last two columns 
of Table [5] we include the observing dates of each spectrum 
used in the present work. 

The observations were made at the 2.15 m telescope 
of the Complejo Astronomico El Leoncito (CASLEO), lo- 
cated in San Juan, Argentina. They were obtained with a 
REOSC spectrograph designed to work between 3500 and 
7500 A and a 1024 x 1024 pixel TEK CCD as detector. The 
spectral resolution ranges from 0.141 to 0.249 A per pixel 
{R = X/SX~ 26400). We refer the reader to CM04 for more 
details on the observations and the data reduction. 

For all the stars, we have several spectra obtained in 
different observing runs. To study the differences in atmo- 
spheric structure with activity level, in this paper we con- 
sider, in most cases, two spectra for each star, chosen be- 
tween those with the better signal to noise ratio. Generally 
we selected the spectra showing the lowest and the highest 
levels of activity, except for HD 177996 and HD 37572, for 
which the least active spectra are very similar to the most 
active ones of HD 22049 and HD 17925 respectively. In this 
way, we built 10 different models. It is important to note, 
given the dependence of activity level with the observation 
time, that all the line profiles used to build the models are 
simultaneous. 

In Figure [T] we show the Scaii index of e Eri obtained 
from our observations (open triangles). The two spectra 
modeled in this paper are represented by full triangles. The 
difference in the Ca II K line fiux between the maximum and 
minimum is 17%. With squares we also present the annual 
average o f the Scaii index. For de tails on the variability of 
e Eri, see lBuccino fc Mauad l|2008l '). 



3 THE CHROMOSPHERIC MODELS 

For each star we built a different chromospheric model, as- 
suming one-dimensional, plane-parallel atmospheres. We si- 
multaneously solved the equations of hydrostatic equilib- 
rium, radiative transfer and statistical equilibrium, using the 
compute r code Pandora. A desc ription of this code can be 
found in lAvrett fc Loesed l|2003l ). 

For a given distribution of temperature with height, we 
self-consistently computed non-LTE populations for 15 lev- 
els of H, 13 of He I, 6 of He II, 15 of Fe I, 8 of Ca I, 5 of Ca 
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Table 1. The stellar sample. Columns 3 to 6 list the stellar parameters (from Ferryman et al. 1997, and from Cincunegui 
& Mauas 2004). The next three columns give the Scall measured by Henry et al. (1996) at CTIO and by Cincunegui &; 
Mauas (2004) at CASLEO, both converted to Mount Wilson Scall compared whit Scall calculated from our models; 
and the last two columns list the observing dates for each spectrum wc used. 



HD (Name) 


S. type 


B-V 


Teff(K) 


[Fc/H] 


■^CTIO 


^CM /•'mod 


^CM /''mod 


Min 


Max 


17925 (V* EP Eri) 


Kl V 


0.86 


4956 


0.10 


0.662 


0.792/0.584 


0.566/0.520 


11/22/02 


12/5/03 


22049 (e Eri) 


K2 V 


0.88 


5110 


-0.14 


0.483 


0.555/0.468 


0.440/O.389 


11/21/02 


3/9/04 


26965 (V* DY Eri) 


Kl V 


0.82 


5203 


-0.25 


0.185 


0.188/0.149 


0.138/0.147 


3/9/04 


8/11/00 


37572 (V* UY Pic) 


KO V 


0.85 


5175 




0.952 


0.703/0.687 


0.691/ — 




11/24/04 


128621 (a Cen B) 


Kl V 


0.90 


5037 


0.24 


0.209 


0.247/0.180 


0.164/0.139 


8/13/00 


9/11/03 


177996 (— ) 


Kl V 


0.86 


5092 




0.861 


0.821/0.798 


0.613/ — 




6/27/02 



II, 7 of Mg I, 6 of Mg II, 21 of Si I, 8 of Na I and 6 of Al I. 
The atomic models we u sed for H and Ca II are described in 
iMauas et all l|l997h and lFalchi fc MauasI (|l998l ). The Ca II 
lines and Lya were computed using Partial Redistribution, 
as it has been done in pre vious chromospheric m od- 
els (like, for example, the lVernazza et al.l (|l98lj ) so- 
lar models) . 

An important element to include in this kind of model- 
ing, in particular for the coolest stars, is the effect of bound- 
bound absorptions due to the numerous atomic and molecu- 
lar lines pr esent in the stellar atm osphere, referred to as line 
blanketing l|Falchi fc Mauaj|l998l ). which plays a crucial role 
in determining both the emergent energy distribution and 
the physical structure of the atmosphere. In solar-type stars 
the most important effects come from neutral or single ion- 
ized metals. In even cooler stars, molecular bands, as CN, 
CO, H2O, etc, could dominate. In this pa per, line blanket- 
ing is treated in non LTE, as explained in lFalchi fc Maua^ 
assuming the source function is given by 

S^=aJ^ + {l-a)B^, (1) 

where is the Planck function and is the mean intensity. 
Q is the s cattering a l bedo, for which we used the expression 
given by lAndersonl l| 19891 ) which depends on wavelength, 
depth and temperature. 

From the finished model, we computed the emitted pro- 
files of H/3 and of the Ca II H and K lines, and modified the 
model until we found a satisfactory match with the observed 
profiles. As a check of the accuracy of the models, we also 
compared the observed computed profiles of the Mg I b and 
the Na I D l ines for each model (de tails of these feat ures can 
be found in lMauas et~aLlll988l and lDi'az et"ai1 l2007'). 

For comparison with synthetic profiles, the observations 
were converted to the stellar surface through 

log{F,^r f/f earth) ^0.35 + 0.4{V + BC)+4log{T,ff), (2) 

where Fsurf is the stellar surface fiux, f earth is the fiux ob- 
served at earth, V is the visual magnitude. BC is the bolo- 
metric correction given byLjohnson ot al. (1966), and T^s is 
the effective temperature for each star, given in Tabled 

Of course, semiempirical models like this one are only 
a first approximation to the structure of stellar chromo- 
spheres, which are neither static nor homogeneous. Regard- 
ing temporal variations, we took care of picking our ob- 
servations at times when no fiares we r e pre sent, using the 
method explained in lCincunegui et al] l|2007l ). Spatial inho- 
mogeneities characteristic of magnetically active stars, like 
starspots or active regions, cannot be resolved on the stellar 



surface. The models presented here, however, can be used 
as a first step t o build two component models as was done, 
for example, bv lMauas fc Falchil l|l996l ). 

Faster temporal variations, like waves, cannot be re- 
produced with this kind of models, of course. We are also 
not considering possible small-scale spatial inhomogeneities 
like, for exam ple, the chrom ospheric bifurcation proposed 
for the Sun bv lAvred (|l98lf ). which should be produced by 
CO cooling. Howeve r, on one hand thi s cooling was prob- 
ably overestimated (|Mauas et al.l Il990l ) , and on the other 
it is probably too slow compared t o atmospheric dynamics 
(|Wedemever-Bohm fc Steffenll2007l ). In any case, homoge- 
neous models provide information on the "mean" state of 
the stellar atmosphere, where the different components are 
weighted by their effect on the emitted radiation, in partic- 
ular on the spectral features under study. 

3.1 Stellar parameters for e Eri 

Before building the model atmosphere, a set of atmospheric 
parameters has to be determined. Both the surface gravity 
and the metallicity are fundamental input parameters in any 
atmospheric model, and the effective temperature, although 
is not needed as input, is used in Ec.[5]to calculate the stellar 
surface fiux needed to analyze the results. 

In Table l3m we summarize several values of these quan- 
tities that can be found i n the literature . Given the astro- 
physical interest on e Eri, iDrake fc Sniitbl (|l993D recognized 
the necessity of determining these parameters with high pre- 
cision and they summarized the methods used to obtain 
them until 1993, and the validity of these determinations. 
To improve these values, they determined the surface grav- 
ity, metallicity and effective temperature in a self-consistent 
way, comparing the equivalent widths of several Fe I, Fe II 
and Ca I lines with theoretical profiles from different model 
atmos pheres. The parameters derived by iDrake fc SmithI 
lll993t) were u s ed in the most recent model for e Eri by 
ISim fc JordanI (|2005l . hereafter SJ05), although they recog- 
nized that the value of log(g) adopted could be too high 
(private communication). 

The difficulty in the calculation of the surface gravity 
is that it is indirectly determined from the values of mass 
and stellar radius. Since these two parameters can be calcu- 
lated more precisely for stars in binary systems, we studied 
another star of our sample, a Centauri B (HD 128621), per- 
taining to the system a Centauri AB. For close systems like 
this visual binary, the stell ar radii and masses can be deri ved 
with an error of 1 to 10% (|Guenther fc Demarauell2000l ). 
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Table 2. Stellar characteristics f or Ep silon Eridani (HD 
22049) from ICavrel de Strobel et alj ll200lh and Table 1 from 
iDrake fc SmithI l|l993^ . 



log(g) 



4.565 

4.61 

4.4 

4.5 

4.5 

4.1 

4.8 

4.19 

4.80 

4.61 

4.75 

4.75 

4.38 

4.57 

4.7 

4.62 



[Fe/H\ 



Reference 



5020 
5000 
5000 
5100 
5040 
5000 
5040 
4990 
5156 
5180 
5000 
5110 
5104 
5135 
5052 



-0.0 
-0.31 
-0.19 



-0.20 
-0.08 
-0.23 
-0.20 
0.05 
-0.09 
0.06 
-0.14 
-0.12 
-0.07 
-0.06 



jirishiia_Swam^ (1966) 
Hearnshaw (1974) 



Oinas (1974) 
Kelch (1978) 
Tomkin fc Lambert (1980) 
Holwcger (1981) 



Steenbock 

^urnashc v (1983) 
Steenbock (1983) 
Abia et al. (1988) 
Bell fc Gustafsson (1989) 
Drake fc Smith (1993) 
Mallik (1998) 
Tomkin fc Lambert (1999) 
Zhao et al. (2002) 
Bodaghcc ct al. (2003) 
AUciidc Pricto ct al. (.2004) 




log m^jg, (g cm"^) 



According tol Cavrel de Strobel et al] (|200ll ). the values 
of log(g) found for a Cen B range from 4.51 to 4.73, with 
an average value of 4.60. We therefore adopted a value of 
log{g) = 4.65 for all the stars in our sa mple, since this 
value is contained in the range given by iDrake fc SmithI 
l|l993f l. considering the error in their calculation {log{g) = 
4.75 ± 0.1). This sa me value of log(g) was adopted by 
iNess fc JordanI l|2008l ) in their recent study of the corona 
and transition region of e Eri. 

Regarding the rest of the stellar param eters, we adopted 
Teff = 5110 K (iTomkin fc Lambertlll999l ). which is close to 
the value bv lDrake fc SmithI ()l993l ) We adopted solar metal- 
licity as a good approximation for e Eri, as has been done in 
all the previous models for this star, since it is a young star 
w hich is probably not me tal defRcient. This was suggested 
bv lKrishna Swainvl ()l966h . who built a grid of model atmo- 
spheres for e Eri with different metallicities to fit the Ca II 
K line and found that using solar metallicity results in the 
best agreement with observati ons. 

In the case of a Cen B, lAvres fc Linskvl (Il976l ) built 
two models for this star assuming in one case solar metal- 
licity and in the other an abundance twice as large. They 
concluded that the computed profiles of the Ca II K line 
differ very little and are in both cases consistent with the 
observations. 



3.2 The model 

To build the atmospheric models for e Eri, as a first step 
we computed a photospheric structure capable of reproduc- 
ing the observed continuum spectrum for this star. Once the 
photospheric model was obtained, we changed the chromo- 
spheric structure to fit the Ca II K and H/3 lines for both 
situations of interest, i.e. the maximum and minimum levels 
of chromospheric activity. This is the first time this sort of 
analysis is made. 

Figure [2] shows the resulting models, which are pre- 
sented in column mass for comparison with the best one- 
component model from SJ05 (their model B). In Figure [3] 
we compare the computed and observed continuum spec- 



Figure 2. Models for e Eri in its minimum (dashed line) and 
its maximum situation (full l ine). For comparison, we show the 
model B from lSim fc JordanI ijiooi) (dotted line) 



trum of e Eri, and in Figure H] and Figure [5] we show the 
comparison of the observed and computed profiles for both 
levels of activity. It is important to note the good agreement 
of the fit, even better than the one by'Thatchcr ct al. (1991!) 
for all the diagnostic lines and continuum. 

On top of the chromosphere, we added a tran- 
sition region with a similar structure to the solar 
one. Since we have no observations of lines formed 
in this region, we could no constrain it further. How- 
ever, the position at which the transition region be- 
gins was adjusted to fit the observed emission of the 
Ca II k Hne. 

There are several differences between our model and the 
one by SJ05. In our model, the photosphere is hotter, the 
temperature minimum region is narrower and the chromo- 
spheric rise has a larger slope. Also, their transition region 
is placed deeper in the atmosphere, i.e., at higher values of 
the column mass. 

The differences between these models may be caused by 
several factors. As we have already noted, SJ05 used a higher 
value of surface gravity which could explain the differences in 
all the thermal structure. The differences in the photospheric 
structure could arise from the fact that we used the complete 
spectra to fit the conti nuum emission, and S J05 used the 
photospheric model by iThatcher" et all (|l99ll ). built to fit 
only the Ca II K line wings, which are formed in the higher 
photosphere. 

Another important factor to consider is the moment of 
the activity cycle in which the observations used to build 
the model were taken. In our case, all the lines used as diag- 
nostics correspond to the same activity level since they were 
all observed simultaneously. But in the model by SJ05, the 
structure of the higher chromosp here and transit i on re gion 
was assembled with the model bv lThatcher et all l|l99ll ) for 
the lower chromosphere and photosphere, whitout taking 
into account that these thermal structures were obtained 
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Figure 3. Comparison between the observed (grey) and com- 
puted continuum (black) for e Eri. 
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Figure 4. Comparison of observed (dashed line) and computed 
profiles (full line) for e Eri in its minimum. 

using line profiles that correspond to different parts of the 
activity cycle. For these reasons, the comparison between 
these models is only qualitative. 

Regarding the differences in the atmospheric structure 
between the maximum and the minimum level of activity, 
the changes occur all along the atmosphere (Figure [2]), from 
the temperature minimum to the transition region. The posi- 
tion of the minimum is the same in both situations, although 
the temperature increases from 3980 K to 4050 K. 

Finally, to check whether our results are affected 
by the adopted value of the metallicity, we computed 
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AHk) 4X(A) 

Figure 5. Comparison of observed (dashed line) and computed 
profiles (full lino) for e Eri in its maximum. 



the emitted profiles for our mode ls with the metal- 
lici ty sivBn by [Z hao et aT] (|2002l ). which was used 
by iNess &: Jordarii (i2008|), and we found no signif- 
icant differences . This result is consist ent with the 
one obtained bv lAvres &: Linskvl (|l976l V 



3.3 The other stars 

To build the models for the other stars in the sample, we 
used solar metallicity and the same surface gravity that was 
used for e Eri. The stellar surface flux was computed with 
Equation[2l using for each star the Tcff values shown in Table 
[2l Since we want to study the changes in thermal structure 
induced by activity, we made the approximation that all the 
stars have the same photosphere than e Eri. 

The models for the less active stars (HD 128621 and HD 
26965 in its maximum and minimum activity level, and e Eri 
in its minimum) are shown in Figure |S] It can be seen that 
all these models have the temperature minimum between 60 
and 100 km higher, and from 20 to 240 K cooler that e Eri 
in its minimum. 

The temperature in the chromosphere, from the 
temperature- minimum region up to 1100 km, increases with 
activity, although the largest differences are in the chromo- 
spheric plateau. These changes with activity are different to 
those obtained for G stars (Paper I) with similar activity 
levels, because in that case only the temperature minimum 
region changed, and the rest of the atmospheric structure 
remained the same. 

An important fact which can be seen in Figure [S] is 
that the differences in the atmospheric structure for a star in 
its maximum and minimum activity levels are comparable 
to the changes seen beetwen two different stars. This fact 
stresses how important it is, when building an atmospheric 
model, the moment at which the observations to be ajusted 
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Figure 6. Models for the less active group. All the models have 
the same structure below 150 km. 
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Figure 7. Models for the more active group. All the models have 
the same structure below 150 km. 



are made, and, in particular, how important it is to use 
simultaneous observations of the diagnostic lines. 

In Figure [13] to [16] we show the observed and computed 
profiles for a Cen B (HD 128621) and HD 26965 in its maxi- 
mum and minimum states. It is important to note the change 
in scale to compare with Figures [J] and [S] since these two 
stars are less active than e Eridani. 

The models for the most active stars (HD 17925 in both 
activity levels, and HD 22049, HD 37572 and HD 177996 
in their maximum) are shown in Figure [7] Again the differ- 
ences in the atmospheric structure for a star in its maximum 
and minimum activity levels are similar to the changes seen 
beetwen two different stars. 

In Figure[7]it can be seen that for the stars in this group 
the temperature minimum is hotter than for the stars in Fig- 
ure[n] and this temperature is almost constant as the activity 
level increases, varying only 50 K. The position of this region 
is also the same for all these stars. The atmospheric struc- 
ture changes with activity everywhere in the chromosphere, 
mainly in the plateau and the rise to the transition region. 

The observed and synthetic profiles for the most active 
stars are compared in Figures[T7]to[2n] Again, it is important 
to note the change in the scale for comparison with the less 
active stars and the good fit in all cases. 



4 NON-RADIATIVE HEATING IN K STARS 



As was mentioned in iKelch et all l| 19791 ), the ratio of the 
temperature in the minimum and the effective temperature 
(Tmin/Teff) givcs au indication of the importance of nonra- 
diative heating in the upper photosphere of stars. In that 
paper, they compare this ratio with Tch to study the trend 
due to spectral type. 

In Figure [S] we plot this ratio versus Scaii, which is 
an indicator of the level of magnetic activity in the chro- 
mosphere for all stars independently of spectral type. The 



Figure 8. Tmin/TLff vs. Scall computed from the models for 

K stars (this paper, squares) and for G stars (Paper I, triangles). 



values of Scaii were obtained by integration of the synthetic 
profiles, and in the figure we include the values obtained 
from the models for K stars built in this paper and those for 
G stars constructed in Paper I. 

In the figure it is possible to observe that there is a sat- 
uration in Tmin- In fact, its value increases with activity up 
to Tmin/Toff ~ 0.79, and after that it remains almost con- 
stant even if activity increases further. On the other hand, 
the computed value of Tmin/Tofr for G stars is larger than 
for K stars with similar activity levels. 

In Figure [9] (left) we show the position of the tempera- 
ture minimum region in column mass as a function of Scaii 
for G (Paper I, triangles) and K stars (squares). For K stars 
the temperature minimum occurs deeper than for G stars. 
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Figure 9. Position (in column mass) of the temperature minimum region (left) and the transition region (right), as a function of 
computed Scall for G (Paper I, triangles) and K stars (squares). 



and there is a tendency for this region to move inward as 
activity grows. In other words, the temperature inversion oc- 
curs deeper for more active stars, indicating that the energy 
deposition starts deeper in the atmosphere as the activity 
level increases, for both spectral types. In Figure |9] we also 
show the position of the transition region (TR), specifically 
the height at which the temperature reaches 36000 K. It can 
be seen that for G stars the chromosphere is more extended 
than for K stars, and that in both cases the TR moves in- 
ward as activity increases. 

To study the energetic requirements to maintain the 
atmospheric structure, we calculated the total net radiative 
loss for each model in the same way as in Paper I. At a given 
depth, the radiative cooling rate <1> (ergs cm~^ sec~^) in a 
given spectral feature (line or continuum) can be computed 
as (Vernazza et al. 1981) 

'i> = 47T J (5, - J,)di/, (3) 

where Si, is the source function and Jv is the mean inten- 
sity at frequency u. A positive value of $ implies a net loss 
of energy (cooling), and a negative value represents a net 
energy absorption. 

Here, we considered line and continua of H, H-, H-ff, 
Mg I and II, Fe I, Si I, Ca II, Na I and CO. The total 
rates for each star are shown in Figure [10] for the less active 
models, and in Figure [TT] for the more active ones. As it is 
expected, the amount of non-radiative energy supplied to the 
chromosphere increases everywhere with magnetic activity. 

In both figures, it is possible to note a region where the 
net cooling rate is negative. This fact was already known for 
the Sun (Vernazza et al. 1981), and was later found in Paper 
I for other G stars, for which negative cooling rates in the 
temperature minimum region were also obtained. Within the 
plane-parallel, homogeneous approximation we are investi- 
gating, this implies either mechanical energy extraction or. 



more likely, that the calculation s have n eglected important 
sources of radiative cooling (see lMauas I l993). 

The main contributions in this zone are H-, Si I, Fe I 
and CO, the same than for G stars. It is important to note 
that since the temperature for K stars is lower in this re- 
gion, there could be an important contribution of several 
molecules which we do not consider in our calculations, like, 
for example, CH, that could act as cooling agents. Consider- 
ing these contributions could bring our computations closer 
to energy balance. 

For the less active models the cooling rate becomes pos- 
itive at around 300 km, implying that there is mechanical 
energy deposition above this height. For the most active 
models, this energy deposition starts deeper in the atmo- 
sphere, i.e. the chromosphere starts deeper. 

Also in the chromosphere, the most important contrib- 
utors to the cooling rate are the same than for G stars, but 
the proportions are difi'erent: for e Eri in its minimum, for 
example, Mg II and Ca II contribute with ~ 9% each, while 
for the Sun these contributions are of 20%. The contribu- 
tion by Fe I, on the other hand, is of ~ 15% in e Eri, but 
only ~ 10% in the Sun. In both cases, almost half of the 
total cooling rate corresponds to line blanketing. 

Finally, to quantify the total amount of mechanical en- 
ergy deposited in the chromosphere, we integrated the net 
radiative cooling rate from the depth in the chromosphere 
where the cooling rate becomes positive to the region where 
the temperature reaches 10* K. To compare the results for 
both spectral types, we normalized the integrated rate, (jiint, 
by the surface luminosity [cT^ff) The resulting quantity, 
therefore, gives an idea of the fraction of the total 
energy emitted by the star that goes into heating the 
chromosphere. The results are shown in Figure [121 where 
it can be seen that there is a unique trend for all stars, in- 
dependently of spectral type. This fact seems to imply that 
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Figure 10. Logarithm of the total coohng rate for the less active 
stars. A positive value of log <& implies a net loss of energy (cool- 
ing), and a negative value represents a net energy absorption. 



Figure 12. Normalized (/lint versus computed Scall index. 
Empty squares represent the K star models from this work and 
full squares indicate the G star models from Paper I. 




O Wgh.(km) 

Figure 11. Logarithm of the total cooling rate for the more ac- 
tive stars. A positive value of log <1> implies a net loss of energy 
(cooling), and a negative value represents a net energy absorption. 



the physical processes that supply the energy to sustain the 
atmospheric structure are independent of spectral type. 

Given the good corelation between Scaii and the nor- 
malized ^int, we fit the data with a polynomial function, 
given by 

= -1-14 lO-'^ + L28 10"* Scaii 
+2.80 10"* sLii - 2.80 10"* Si^,„ . (4) 

In Figure [12] it can be seen that the fit is very good and, 
therefore, the energetic requirements of a given star can be 
estimated from its chromospheric activity level as mesured 
by Scaii. 



5 DISCUSSION 

One of the main goals of chromospheric modelling is to ac- 
curately estimate the radiative losses in the chromosphere 
in detail, using only the information that can be obtained 
from the observations, without any assumption about the 
physical processes involved. In this way, these losses can be 
equaled to the energy requirements that any proposed mech- 
anism of chromospheric heating should match. 

For example, we saw in the previous section that 
the contribution of the different features to the total 
cooling rate is not the same for G and K stars. In 
particular, the Ca II, Mg II and Fe I relative con- 
tributions are not the same for both spectral types. 
Therefore, it might not be correct to scale the rela- 
tive contributions computed for the Su n to K stars, 
as it has been done somet imes (see ICuntz et al.l 
I1999I . and iRammacher et aD[2005l '). 

Cuntz et al. (1999) computed theoretical two- 
component models for K dwarfs of different activity levels. 
They proposed that the energy is deposited in the chro- 
mosphere by acoustic and magnetic shocks, and found that 
these shocks are stronger and are produced deeper in the 
chromosphere as the activity of the star increases. This re- 
sult is in agreement with our calculations, which shows that 
the energy deposition is larger and deposited deeper with 
increasing activity. 

On the other hand, they reproduced the lineal trend 
between the Ca II H and K lines fiuxes and the rotational 
period, although their computed fluxes are smaller than the 
observations, which could be due to their sketchy calculation 
of the radiative cooling rate. 



6 SUMMARY 

In this paper we present chromospheric models for six K 
dwarfs, including e Eridani, with similar photospheric prop- 
erties but different magnetic activity levels. In most cases 
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we computed models for two moments of the activity cycle 
for the same star. 

These models were based on, and reproduced very well, 
the Ca II H and K and the H/3 line profiles for all the stars in 
our sample. The reliability of the stellar atmospheric models 
was checked with other features, the Na I D and Mg I b lines. 
Also for these lines we found very good agreement between 
computed and observed profiles. 

We found that the changes in atmospheric structure in 
K dwarfs with activity arc produced all along the chromo- 
sphere, from the region of the temperature minimum to the 
transition region and mainly in the chromospheric plateau, 
independently of the activity level of the star. This was not 
the case for the G dwarfs modelled in Paper I, since for the 
less active G stars the changes with activity occur only in 
the region of the temperature minimum. 

The ratio of the minimum and effective temperatures 
(Tmin/Toff) can give an idea of the importance of non- 
radiative heating in the upper photosphere of stars. Both 
for K and G stars, this value increases with activity up to 
Tmin/T'eff ~ 0.79, wlicrc it saturates, and it remains con- 
stant even if the activity level increases further. On the other 
hand, the computed value of Tmi-a/Tgs for G stars is larger 
than for K stars with similar activity. 

For both spectral types, the position of the temperature 
minimum moves inward as activity increases, implying that 
the chromosphere starts deeper in more active stars. This, in 
turn, implies that as the activity level increases, the energy 
deposition occurs deeper in the atmosphere. 

On the other hand, the transition region is placed at 
higher column masses for G stars than for K stars, and in 
both cases it moves inward as activity increases. 

Regarding the energetic requirements, the integrated 
chromospheric radiative losses, normalized to the surface 
luminosity, show a unique trend for G and K dwarfs when 
plotted against Scaii, the main proxy of stellar activity. This 
might indicate that the same physical processes are heating 
the stellar chromospheres in both cases. We calculated an 
empirical relationship between the Scaii index and the en- 
ergy deposited in the chromosphere, which can be used to 
estimate the energetic requirements of a given star knowing 
its chromospheric activity level. 

There axe significant differences in the contributions of 
Mg II, Ca II and Fe I to the total net cooling rate in the 
chromosphere between G and K stars, which implies that 
values obtained for a given star should not be extrapolated 
to another one of a different spectral type. In both cases 
about half of the total rate is due to line blanketing. 
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Figure 13. Comparison of observed (dashed line) and computed 
profiles (full line) for a Can B (HD 128621) in its maximum. 
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Figure 14. Comparison of observed (dashed line) and computed 
profiles (full line) for a Cen B (HD 128621) in its minimum. 
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Figure 20. Comparison of observed (dashed line) and computed 
profiles (full line) for HD 177996. 
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